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Abstract 

We propose a novel protocol for storage and retrieval of photon wave packets in a A-type atomic 
medium. This protocol derives from spectral hole burning and takes advantages of the specific 
properties of solid state systems at low temperature, such as rare earth ion doped crystals. The 
signal pulse is tuned to the center of the hole that has been burnt previously within the inhomo- 
geneously broadened absorption band. The group velocity is strongly reduced, being proportional 
to the hole width. This way the optically carried information and energy is carried over to the off- 
resonance optical dipoles. Storage and retrieval are performed by conversion to and from ground 
state Raman coherence by using brief vr-pulses. The protocol exhibits some resemblance with the 
well known electromagnetically induced transparency process. It also presents distinctive features 
such as the absence of coupling beam. In this paper we detail the various steps of the protocol, 
summarize the critical parameters and theoretically examine the recovery efficiency. 
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I. INTRODUCTION 



Future long distance quantum information networks are expected to involve quantum 
repeaters. To this end, intense investigation is going on all over the world to interface a 
quantum state of light with an atomic ensemble. In protocols based on electromagneti- 
cally induced transparency (EIT), the input signal is adiabatically converted into an atomic 
ground state coherence, without excitation of optical coherences. This class of protocols 
has been investigated intensively both theoretically [ij^, 3] and experimentally, leading to 



storage and retrieval demonstration of both discrete 4, ISl] and continuous [g, [3] quantum 
variables. As it bypasses the optical coherence excitation, the EIT based protocol is well 
adapted to atomic vapors, with strong transitions and short optical dipoles lifetime. 

Rare earth ion doped crystals (REIC) also offer long ground state coherence lifetime, 
which makes them attractive for quantum storage. However, in contrast with the atomic 
vapors considered in EIT experiments, they also exhibit long optical coherence times. In 
addition, the active ions are strictly motionless, the hyperfine sublevel relaxation time may 
reach tens of hours, and the optical transition frequency is distributed over a very broad 
inhomogeneous width. 

Because of these specific features EIT-based protocols may not be the most appropriate 
approach to quantum storage in this type of medium. Actually, different protocols have been 
considered recently, based on rephasing processes 

,[io|,[n[[i2|,Q- 

In the present paper. 



we propose a new adiabatic protocol, taking advantage of the specific REIC properties. 

It is well known that a light pulse slows down when it propagates through a transparency 
window (TW) in an absorbing medium. Deeper and narrower the transparency window, 
slower the light pulse. The group velocity v is of the order Ao/ao, where Aq and ao re- 
spectively stand for the TW width and the linear absorption coefficient outside the TW. 
Provided the input signal spectrum is narrower than the transparency window, the pulse 
undergoes little energy loss and shape distortion. 

In EIT, the TW is opened by a strong classical coupling field, interacting with the atoms 
simultaneously with the input signal. As the input signal is slowed down to group velocity 
V « c, the pulse length and energy are reduced by a factor v/c « 1. The pulse energy 
is carried over to the coupling field, and the input signal quantum state is mapped into the 
ground state Raman coherence. This coherence behaves as a spin wave that adiabatically 
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follows [l^ the signal pulse, at the same velocity v. By switching off the coupling beam 
one closes the TW. This simultaneously destroys the slowly propagating pulse and freezes 
the spin wave that already contains most of the signal information 15|. To retrieve the 
stored information, one turns on the coupling beam. Interacting with the coupling beam, 
the frozen spin wave restores the signal pulse. Conversely, the reviving light pulse drives 
the spin wave to the medium output where the initial energy is finally returned to the light 
signal by the coupling beam. 

In REIC, relying on the homogeneous width narrowness, the hyperfine relaxation slow- 
ness, and the absence of translational motion, one can prepare a narrow TW, by frequency 
selective optical pumping to a ground state shelving sublevel. This is nothing but the well 
known spectral hole-burning process. With an hyperfine level lifetime ranging from a few 
seconds to several hours, the preparation step can be completed long before the arrival of 
the input signal. After the preparation step, an input signal that is sent through the spectral 
hole undergoes the same slowing down as in EIT 16|]. In contrast with EIT, the energy and 
the signal state are both carried over to the optical coherence of the off-resonance atoms. 
The transient optical dipoles are created along the path of the signal field. Their natural 
lifetime is much longer than the pulse duration, but they vanish adiabatically as the light 
pulse passes by. To "stop" the light and store the information, we propose to convert those 
optical dipoles into Raman coherences, while the slowed down pulse is still confined within 
the absorbing medium. This can be performed by an auxiliary vr-pulse, much shorter than 
the inverse TW width. Another vr-pulse shall restore the optical dipoles and release the 
stored information. 

The paper is arranged as follows. We describe slowing down and storage in section II. 
Section III is devoted to the signal recovery. The latter section represents the core of the 
paper. We detail the revival of the signal from the restored optical coherence. Specifically 
we determine the rising time of the retrieved signal, and its temporal profile. We finally 
compute the recovery efficiency, taking account of the size of the absorbing medium and of 
the signal spectral width. 
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FIG. 1: (color online) inhomogeneous distribution G(A)/Go with a Ao-wide spectral hole. The 
inset represents the A-shape three-level system. The input signal pulse excites transition a-b. 

II. SLOWING DOWN AND STORAGE 
A. The slowing down step 

We consider an ensemble of motion-less three-level atoms in a A configuration. The 
excited state \b) is connected by optical transitions to the ground substates \a) and |c). The 
light pulse to be stored interacts with transition \a) — \b) (see Fig. [1]). The inhomogeneous 
broadening of the optical line is much larger than the pulse width. A hole is burnt in the 
absorption profile over an interval Aq around frequency ujq. This is achieved by pumping 
the corresponding atoms to some additional long lifetime shelving level. The homogeneous 
linewith 'jat is assumed to be much smaller than Aq. All these assumptions are consistent 
with the properties of some rare earth ion doped crystals at low temperature. The absorbing 
atoms are distributed as a function of the transition frequency Uab- Since 'jab << Aq various 
shapes can be given to the burnt hole. 

Slow light with persistent hole burning was already discussed thoroughly [l^ and we just 
summarize the main features relevant to the proposed storage protocol. The incoming field 
E{z,t), centered at frequency c^o, propagates along direction z and is expressed in terms of 
the slowly varying envelope A{z,t) as: 




(1) 



The pulse propagation is described by the linearized wave equation: 



d,Aiz,t) + -dtAiz,t) 



i—V{z,t) 



(2) 



The polarization density is expressed as the sum V{z,t) = fiat J dAG{A)aab{A; z,t) of the 
off-diagonal density matrix element 



oo 



aab{A; z, t) = -i^ drAiz, t - r)e(^^-^-)^ (3) 

where fiab stands for the transition dipole moment. The sum is performed over the distribu- 
tion G{A) of the atomic transition frequency Uab = ujq + A. The distribution G{A), initially 
uniform over a broad interval, has been depressed by hole burning over a Ao-wide spectral 
window centered at A = 0, where G{0) = 0. 

In the weak field limit, the equation is most conveniently solved in spectral domain in 
terms of linear susceptibility. Time-to-frequency Fourier transform of Eq. [2] leads to: 

d.Aiz, n) + i-Aiz, n) = -i—V{z, ^) (4) 
c eo 

where: 

V{z,n) = ]^e,x{^, + n)A{z,n) (5) 

and 

X(.„ + «)^|^/.A^|^ (6) 

heQ J A + i-iab 

The absorption coefficient can be expressed in terms of the imaginary part of the 

susceptibility as a{uj) = —x"{i^)/k. At large spectral distance from the burnt hole, this 
quantity reduces to: 

«o = vrGo^, (7) 

where Go stands for the uniform inhomogeneous distribution away from the hole. 

To get more tractable expressions in the burnt hole region, we expand the susceptibility 
to second order in Q. Then Eq. H] reduces to: 



d,A{z, n) + 



I— + -Oi{uJo + il) 



Aiz,n) = (8) 



where the absorption coefficient and the inverse group velocity are respectively expressed 
as: 

a(L,„ + !J) = a„ (jdAg{A)C(A) + ^ J dA?^^aA)\ (9) 
i4.^/.A,(A)^^^ (10) 
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where g{A) = G{A)/Go denotes the normahzed inhomogeneous distribution and C{A) = 
7r"^7af,/(A^ + lib) represents the Lorentzian homogeneous hne profile. 

If the incoming pulse spectrum is much narrower than Aq, absorption is well described 
by the f2-independent coefficient ci;(ti;o)- Then, by Fourier transforming Eq. [H] solution back 
to time domain one obtains: 

A{z, t) = A{0, t - ^/t;)e-^°("«)" (11) 

This is a travelling wave propagating without distortion at velocity v through the L-long 
slab. According to Eqs. [9] and \T0\ a^Uo) and 1/v are respectively of order ao7afe/Ao and 
ao/Ao. The wave is not attenuated provided ao-^7ab/Ao << 1. With respect to free space, 
the wave envelope is spatially compressed by factor v/c when light travels inside the atomic 
medium, the pulse duration T remaining unchanged. 

The entire pulse may be confined inside the medium if the group delay L/v is larger than 
T, which also reads AqT << aoL. This condition must be consistent with the width of the 
burnt hole, since a lower limit to the pulse spectral width is given by 1/T. To account for the 
finite size of the transparency window we use the second order term in Eq. [9l At the hole 
center, d\g{A) ^ Aq ^. Then, with 'jab « Aq, the transmission factor at spectral distance 
Q from the hole center reads as exp{—^aQLQ'^ / Aq) . Hence, the pulse spectral width is much 



smaller than the transmission window width if AqT >> yoo^- The two conditions can be 
combined into: 

(12) 

showing that the optical density is the only critical parameter. Not surprisingly, a similar 
condition has been derived for the EIT process 

To check the validity of the susceptibility second order expansion, let us consider a burnt 
hole with a gaussian profile. Then, provided ■jab << Aq, the susceptibility reads as: 



X{uJo + i2) = — 



1-e ^o]+'f{^] 
Ao 



(13) 

where F(x) represents the Dawson integral. To second order in Q, Eq. [13] can be expanded 



as: 



xi^o + n) = -{——-t-^ (14) 



k \^/^T Ao A 



As illustrated in Fig. [21 the propagation of a gaussian pulse, as computed to second order 
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FIG. 2: propagation of a gaussian pulse through a gaussian spectral hole. The bold, solid and 
dashed lines respectively represent the initial pulse shape, the exactly calculated transmitted pulse 
shape and a second order susceptibility approximation of the transmitted pulse profile. Opacity is 
aoL = 100 and the displayed profiles correspond to AqT = 5 and AqT = 10. 

approximation with the help of Eqs. [Eland [HI keeps very close to the exact solution derived 
from Eqs [H and [131 



B. The storage step 

Being spatially compressed in the medium, the signal pulse looses most of its energy. If 
the medium is comprised of the active atoms alone, the dielectric host being ignored, the 
index of refraction is unity at uq and the field amplitude is conserved at the input and 
output sides of the slab. The spatial density of energy is also conserved. By summing this 
density over the compressed envelope, one can easily show that the optically carried energy 
is reduced by the factor v/c. The same phenomenon occurs in EIT. In this process energy 
flows from the incoming pulse to the coupling beam at the input side of the medium and is 
restored by the coupling beam at the output side. In the hole burning process the situation 



ho 



is different since there is no coupling beam. Then the atoms located outside the burnt 
serve as a temporary energy tank. As discussed a long time ago by D. Grischkowsky 
this transient off-resonance excitation can be analyzed in terms of adiabatic following. 
With the condition limf^_oo A{z, t) = 0, repeated integration by parts of Eq. [pleads to: 

1 ,, °o / _ • \ "+1 an 
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FIG. 3: (color online) time diagram for storage and retrieval. Incoming and restored signal temporal 
profiles are represented in the input and ouput sides of the memory medium. The incoming pulse 
is tuned to the burnt hole center, interacting with the atoms on transition \a) — \b). The storage 
vr-pulse is applied at time ^,^1 = L/{2v) on transition |6) — |c), when the incoming pulse is halfway 
from the slab ends. Retrieval is accomplished by a second vr-pulse at time t,r2- The restored profile 
is represented in red. In the absence of storage pulse, the original pulse emerges from the slab with 
the delay L/v (magenta line). The same color code is used in Fig. [5l 



To lowest order in fiab-^iz,t)/h/A, off-resonance atom coordinates in the Bloch sphere equa- 
tor read as: 



u{A;z,t) 



.v{A;z,t) 



(16) 



h A n A2 

where (Jafe(A; z, t) = [n(A; z, t) + iv{A] z, t)] /2. Therefore, despite the weakness of the driv- 
ing field, the off resonance optical dipoles with different detunings do not dephase from each 
other, as long as the field is present. On the contrary, in the absence of the field, the position 
in the Bloch sphere equator would rotate at angular velocity A. 

The time diagram of the proposed storage and retrieval protocol is illustrated in Fig. [31 
By shining a brief vr-pulse on transition b-c one can convert the optical dipoles into Raman 
coherence. Illumination by the vr-pulse must take place at a moment tjri when the signal 
pulse is confined in the medium. A pulse of duration Tc successfully converts dipoles up to 
1/tc from resonance. This way one decouples the atoms from the signal pulse. Since the 
compressed signal pulse carries a small part of the initial energy, the information is really 
frozen inside the medium. Deprived of the co-propagating dipole wave, the signal pulse 



8 



recovers the velocity of light in vacuum and exits the medium. Indeed it takes a time of 
order I/Aq, much larger than L/c, to rebuild the macroscopic polarization needed to slow 
down a light pulse. The pulse accelerates to c as a whole, which means its duration is 
reduced by the factor v/c and its spectral width is conversely magnified by the factor c/v. 



III. RECOVERY OF THE STORED SIGNAL 
A. Expression of the retrieved signal 

In order to restore the signal field one applies a second brief vr-pulse at time tT^2 (see Fig. 
|3]). Calculating the restored signal represents the core of this paper. As discussed in section 
II, the atomic dipoles are locked to the incoming pulse in the first stage of the process. When 
they are restored by the second 7r-pulse, they come to life in a field-free medium. Since no 
field is yet able to lock them, they rapidly dephase from each other. At the same time, they 
give rise to a radiative response that will be able, ultimately, to keep them in phase. One 
may wonder whether dephasing may affect the restored signal amplitude. 

At the moment, let us assume that the conversion operated by the two vr-pulses is perfect. 
The recovered polarisation can be expressed as: 

Pi(^, t) = -i^Yit - K^) j dAG{A) drA^z, - r)e(^^-7-)(-+*-*-) (17) 

where Ain{z, t) represents the initial field amplitude and Y{t) denotes the unit step function. 
The polarisation also includes the atomic reaction to the restored field Aout{z,t): 

V2{z,t) = -i^ j dAGiA) J^^ drAoutiz,t-T)e^'^-''^^^^ (18) 

The wave equation for Aout{z,t) reads as: 

1 k 

d,Aout{z, t) + -dtAoutiz, t) = -%- [Pi(2;, t) + V2{z, t)\ (19) 
c eo 

In the same way as in EqslHEl Fourier conjugation leads to: 

/r2 1 \ ~ k ~ 

d.Aoutiz, + -eox(^o + ^) Aout{z, n) = -i-Vi{z, fi) (20) 

V c / / eo 

Solving the equation and transforming back to the time domain, one obtains: 

Aoutiz, t) = -^J dAg{A) jf dr rfr'e('^"^-)(-+-')7^(t.l -r^t-t^^- r', z) (21) 
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where 



7^(^^ t, z) = dz'An{z\ t') j ^e*nt-<§+ifcx(-o+n))(.-.') (32) 



B. Restored signal rising time 

Some time is needed for the signal to revive. First we verify that it starts from zero at 
t = t,r2- We specify the hole profile assuming g{A) = 1 — exp(— A^/Aq). As illustrated in 
Fig. [21 the susceptibility can be reduced to its second order expansion in Q. Then, as shown 
in the Appendix, the restored signal can be written as: 

AUz,t) = [ dAg{A) r dr dr'e(^^-^-)(^+^') 

(23) 



T- 

n=0 



where 6 = t — — t' ■ This expression vanishes at t = tT^2 since the interval of integration 
over r' reduces to 0. Next we determine the signal rising time. 

The restored signal is observed at the slab output, at z = L. The n = term dominates 
in Eq. [23] provided: 

An(2;',t.i - r) » /3^dl [Ainiz',Ui - t){L - Z%,^^_^g (24) 

^0 

Let Ain{z, t) be a Fourier transform limited pulse of duration T. Given that the pulse spec- 
trum is assumed to be narrower than the burnt hole, and that d'^Ain{z, t) ~ Ain{z, t) / (t'T)^, 
the condition reduces to: 

Ao min {t - U2, L/v) « {AoTf (25) 
Within the time interval when this condition is satisfied, the restored field is given by: 

dAgiA) / dr / rfr'e(*^-^-)(^+^'U,„(L - vO, - r) 

Jo Jma,x{0,t-t^2-L/v) 

(26) 

The sum over r' indicates that only slices at distance < v{t — 1(7^2) contribute to the restored 
signal. Let us first focus on the time region t < tT^2 + L/v. Then the lower boundary of the 
integral over r' is 0. Significant contributions only come from r and r' values smaller than 
Ap ^ Given that AqT » 1, we can consider that Ain{L — v9, t-^i — r) does not vary on this 
timescale. Then, taking Ain{L — vO^It^i — r) out of the integrals at r = r' = one obtains: 

Aout{L, t) = An{L - v{t - U2),Ui)k [{t - t.2)Ao/2] (27) 
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FIG. 4: Evolution of the retrieved signal amplitude after illumination by the second vr-pulse at 
t^2- {^) infinite bandwith operation, (b) finite bandwith operation: Ai = 5Aq. The wavy line 
appearance reflects a sharp cut-off in the frequency domain 

where 

n{x) = / dAg{A) / dr / dr'e^''^-^'^''^^^^^'^ (28) 

27r J Jo Jo 

= (^1 - ^) (l - e-"' + xy^ [1 - erf (x)]) (29) 

As illustrated in FigJH the field recovery is 90% completed at time t = tn2 + 2/Ao. 

Thus far we have assumed that storage and retrieval are accomplished with an infinite 
bandwidth, all the atoms contributing to the process, whatever their distance A from the 
burnt hole. The atoms contribute with a weight 1/A^ to the restored signal, as can be seen 
by considering the (t — t .,,2)^.0 >> 1 limit of Eq. [271 

Aout{L,t) = A,n{L - v{t - t^2),t.i)^ / dAg{A) (30) 

27r J [A + i-fabr 

Therefore, finite bandwidth operation should reduce the recovery efficiency by a factor close 
to 1 — Ao/Ai, where 2Ai denotes the process bandwidth. With the same gaussian hole as 
above, this factor equals erf (y) — -^[1 — exp(— y^)]/?/ where y = Ai/Aq. The corresponding 
time evolution of the restored field is shown in Fig. H] for Ai = 5Aq. The wavy line 
appearance results from the sharp cut-off in the frequency domain. 

C. Established signal 

Although the series expansion in Eq. [23] was convenient for studying the signal birth, 
this expression is difficult to handle when Eq. [25] is no longer satisfied. A different approach 
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is then needed. Once the signal is estabhshed, the quantity TZit-^i — T,t — t,,-2 — t' , z) evolves 
on a timescale much longer than Aq ^ and can be taken out of the integral at r = r' = in 
Eq. [211 The restored signal amplitude reduces to: 

Aout{L, t) = (l- ^) n{Uu t - U2, L) (31) 

To go further we have to specify the incoming pulse profile. We proceed with the gaussian 
field: 

An(0,t)=e-^ (32) 
According to Eq. [HI propagation from the input side to z' changes the field amplitude into: 



AnT 



t- 



A^n{z',t) = ° =e ^(^o^'+"o^') (33) 

Substituting Eq. [2S1 hito Eq. [22] and using the susceptibility second order expansion one 
obtains: 

^Ao'Ao'Ao^ Jo ^2Txa{p~u)^/SlT^ + au ' 
where a = ^pn . To optimize recovery, we assume that storage occurs when the incoming 
pulse sits halfway from the slab ends, which means t^^x = L/(2v). In order to satisfy the 



double condition v^ooX << AqT << aoL, which is needed to confine the pulse within the 
slab both spatially and spectrally, we set AqT = b(aQL)^^^, where b is close to unity. In 
Fig. [5] we display the time evolution of both the original and the restored signals for various 
values of the optical density. The upper boxes represent the restored field as a function of 
t — tvr2 at the slab output. The lower boxes represent the initial field as a function of t — t^ri, 
in the absence of storage and recovery pulses. The computation is performed with b = 0.6. 
In the figure we indicate the size of T for each value of aoL. Propagation through the 
burnt hole affects the pulse spectral wings, stretching the temporal profile. The size of the 
corresponding elongated pulse width, denoted Tg and given by = T [1 + a^L/ (AqT)^] 1/2, 
is also indicated in the figure. As expected, as a^L is increased, a larger part of the incoming 
signal is comprised within the temporal storage window L/v. In other words, as oqL is made 
bigger, a smaller part of the initial field leaves the slab before the storage time t^^i, or has 
not yet entered the slab at time t^^i. With growing oqL, the restored field is less distorted 
with respect to the incoming signal. Finally, the restored field drops rapidly at time longer 
than L/v, since the finite size of the slab truncates the trailing edge of the incoming pulse. 
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FIG. 5: (color online) original (magenta lines) and restored (red lines) field as a function oi t — t-^i 
and t — respectively for different values of aQL. The original pulse is gaussian and reaches 1, 
its maximum value, at the origin of time at the slab input. Storage is performed at tyri = L/{2v). 
One adjusts the pulse duration to the opacity of the medium by setting: AqT = 0.6(ao-^)^^^- The 
original and stretched pulse durations T and Tg are indicated for each a^L value. The group delay 
L/v determines the temporal memory depth of the slab. 

It should be noticed that Eq. [26] is not valid only in the leading edge of the restored 
signal. According to Eq. [25l this simple form of the retrieved field is valid for any value of 
i — tn2 provided AqL/v « (AqT)^. Since, according to Eq. [221 the signal mainly comes 
from depth v{t — ^,^2), it rapidly vanishes when v{t — ^,^2) exceeds the slab dimension L. 



D. Recovery efficiency 

The recovery efficiency is simply defined as the ratio of the incoming and restored energies: 

JZdt\AUL,t)\' 

This quantity is displayed in Fig. [H as a function of \JolqL. The computation is performed 
with the same parameters as in Fig. [5l The expected recovery efficiency behavior is similar 
to what is observed in EIT. 

The question of the control field dynamic adjustment in EIT experiments was raised 
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FIG. 6: recovery efficiency as a function of yJa^L. The original pulse is gaussian. Storage is 
performed at tj^i = L/{2v). One adjusts the pulse duration to the opacity of the medium by 
setting: AqT = 0.6(ao-^^)^/''- 



recently [17[. It was shown that the temporal variation of the control field at writing and 
retrieval could be optimized to improve the recovery efficiency. According to these works, one 
can map an arbitrary input pulse on the optimal spin wave, unique for each opacity value, 
by appropriate temporal shaping of the control field during the writing step. Conversely, the 
stored information can be retrieved in any desired temporal shape with adequate profiling of 
the readout control field. In the present situation we apparently lose this degree of freedom 
since the burnt hole cannot be reshaped during the storage process. However, the z spatial 
coordinate plays a role similar to the time coordinate. A z-dependent spectral hole profile 



might act on the input pulse in a similar way as the time dependent control field in 



Such a ^-dependent profile can be achieved by transverse illumination of the slab 



T. 



IV. CONCLUSION 



Relying on the specific properties of REIC, we have proposed a new optical storage 
protocol that is based on persistent spectral hole burning. As in EIT, the input pulse 
propagates with negligible absorption and is adiabatically coupled to a collective atomic 
state coherence. In contrast with EIT, inhomogeneous broadening plays a crucial role and 
information is initially carried over to off-resonance optical dipoles. Storage relies on the 
conversion of optical dipoles into Raman coherence and back with the help of brief auxilliary 
TT-pulses. Requirements are similar to those of EIT, and the expected efficiency is also 
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similar. With respect to EIT, the main practical feature is the absence of coupling beam. 
The restored signal is temporally separated from the read out 7r-pulse. Another interesting 
feature is the spectral separation of the atoms that carry the information and those that 
are involved in the hole burning process. We have derived a semi-classical theory but the 
extension to quantum storage would be straightforward. 

Acknowledgement 1 The authors are pleased to acknowledge stimulating discussions with 
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APPENDIX A: RESTORED FIELD CALCULATION 



This calculation aims at using the properties of the Dirac peak and its derivatives. With 
an inhomogeneous distribution given by g{A) = 1 — exp(— A^/Aq), the sum over Q in Eq. 
can be expressed in terms of the susceptibility second order expansion as: 



2tt J 271 ^ ■ ^ 

where = (l — f)- To extract the Dirac peak one expands the exponential in the 
following way: 

dn 



dn in(^t~^)~t3^^ ^ 1 

271 



E 

n=0 



' (3 z-z'' 
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With the help of the Dirac peak derivative definition: 



dn 



rf^" r dn 



27r^ ^ (it2n / 271 



and of the property: 



2n 



(A.3) 
(A.4) 



—AUZ', t') ( \ 5(2") ( t - ^ - ^^^^2n 

V \ V / \ V 



one can perform the sum over z' in Eq. [22] as: 

v-iio,.ut)dinAM,t'){z-zr].,=...t 

(A.5) 

where Ilx,y{t) = 1 when t G [x,y] and vanishes outside. Finally, substituting Eq. IA.5I into 
Eq. [22] one obtains: 



•Aouti^^i t) 



' 27X 



dAg{A) / dr 



max(0,t— t7r2— ^Z") 



oo ^ 
n=0 



(A.6) 



dl?[An{z\t^i-r){z-z'rU 
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where 6 == t — t.„2 — t' ■ 
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